The dynamic thermal expansivity of liquids near the glass transition. 
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Based on previous works on polymers by Bauer et al. [Phys, Rev. E (2000)], this paper describes a 
capacitative method for measuring the dynamical expansion coefficient of a viscous liquid. Data are 
presented for the glass-forming liquid tetramethyl tetraphenyl trisiloxane (DC704) in the ultraviscous 
regime. Compared to the method of Bauer et al. the dynamical range has been extended by mak- 
ing time-domain experiments and by making very small and fast temperature steps. The modelling 
of the experiment presented in this paper includes the situation where the capacitor is not full be- 
cause the liquid contracts when cooling from room temperature down to around the glass-transition 
temperature, which is relevant when measuring on a molecular liquid rather than polymer. 



The glass transition occurs when the configurational degrees of freedom of a liquid are frozen in. Below the glass 
\ transition temperature, Tg, only isostructural contraction takes place as temperature is decreased further. The mea- 

00 ' sured thermal expansion coefficient, a (and heat capacity, c„) are therefore lower in the glass than in the equilibrium 

liquid. This change of the thermal expansion coefficient (and the heat capacity), is probably the most classical signa- 
ture of the glass transition, and a figure illustrating this change (see Fig. [1} is almost inevitably the starting point of 
introductory talks or texts on the glass transition (eg. Ref. [lQ)- 

The change in the heat capacity at the glass transition Ac p = c p/ ;/ 1? — Cp r gi ass has been studied extensively and is 
widely believed to play a role for the dynamics of liquids close to the glass transition. The change in expansion 
coefficient Aa p = Cipjjq — CCp.giass nas received less attention, but is of similar importance. This is seen, for instance, in 
the literature related to the Prigogine-Defay ratio, a dimensionless number characterizing the glass transition 3 - - — . 
The glass is an out-of -equilibrium state and therefore the values of the thermodynamic derivatives are not rigor- 
' ously well defined. They depend on cooling rate and also on the time spent in the glassy state. Contrary to this 
the linear response of the metastable equilibrium liquid state is well defined and history independent 7 . The linear 
expansion coefficient of a viscous liquid close to its glass transition is dynamic, that is time (or frequency) dependent 

1 ^ i , with short times giving a low (glass-like) value, 0Cpj as i, while long times give a higher liquid value, 0Cp iS i ow . The differ- 

ence between these two levels Aa.pjj n = 0Lp iS \ ow — Oipj as t thus gives well-defined information on the configurational 
Tj" | part of the expansion coefficient. Likewise Ac p /,„ = Cp iS \ 0W — Cpj ast is well defined. 

The relaxation between the fast and the slow response takes place on a certain time scale which is temperature 
dependent. Considered in this way the measurement of the expansion coefficient can be viewed as a type of spec- 
troscopy, which gives both a relaxation time and a spectral shape analogous to other methods like dielectric spec- 
troscopy or mechanical spectroscopy. The study of the temperature dependence of relaxation times and of the spec- 
""• , tral shape of different response functions is vital for understanding the viscous slowing down. There is a general 

■ belief that the liquid has a relaxation time, which is fairly well defined independent of probe, but also suggestions 
, that different processes may decouple from each other at low temperatures 8 . 

■ There are good scientific reasons to study the dynamic linear expansion coefficient, but almost no data of this 
L" , type are to be found in literature. The time-dependent expansion coefficient can be found by studying the change in 

. in ' volume as a function of time after a temperature step. Such volume relaxation experiments are very classic in glass 
^ \ science and still important^ - —. However, volume relaxation experiments are traditionally performed as non-linear 
5-j ■ aging experiments, i.e., with large amplitudes in the temperature jump. This type of experiment gives information 
- y.' on the relaxation of the configurational degrees of freedom, but the expansion coefficient and its characteristic time 
scale cannot be determined, because the results depend on the amplitude and sign of the temperature jump. For 
sufficiently small temperature steps this is not the case, defining the linear response regime. 

The only linear dynamic data we are aware of were reported about a decade ago by Bauer et alM^ followed 
by a paper by Fukao and Miyamoto 15 . These papers reported frequency-domain measurements on thin polymer 
films, performed with temperature scans at a couple of fixed frequencies, covering 1.5 decade of the dynamics. 
The measurements were pioneering, but 1.5 decade is not very much for studying relaxation in viscous liquids, 
because the relaxation is extremely temperature dependent and quite "stretched", which means that even at one 
fixed temperature the relaxation covers several decades. 

The technique developed by Bauer et al. is based on a principle where the sample is placed in a parallel plate 
capacitor such that it is the sample that maintains the spacing between the plates. Thus a change in sample volume 
in response to temperature change leads to a change of the capacitance. This principle is also used in the present 
work. The advantage of this technique is that capacitance can be measured with high accuracy and it is this accuracy 
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which makes linear experiments possible. 

The use of sample-filled capacitors for measuring an expansion coefficient is not unique and it has been done 
by others before and after Bauer et al. (see eg. Ref. Il6l - [l9|) in capacitative scanning dilatometry, i.e. working in a 
temperature ramping mode. Capacitative scanning dilatometry has to our knowledge never been used on simple 
liquids. It is particularly useful for studying thin polymer films because the signal gets better with a thin sample. 
The technique has been used for determining the glass-transition temperature for example as a function of film 
thicknes o 17 ' 18 or as a function of cooling rate 1 -. The main focus of these papers is on the temperature dependence 
of the expansion coefficient, while little attention is given to the absolute values. There have been no studies of the 
dynamics since the pioneering work of Bauer and no attempts to extend the dynamical range. 

To the best of our knowledge there are no measurements of the dynamic linear expansion coefficient of molecular 
liquids. The reported data from scanning dilatometry and non-linear volume relaxation are also mainly for polymers, 
while data on molecular liquids is relatively scarce. This may be due to the higher technological importance of 
polymers. It is probably also related to the fact that working with molecular liquids requires other experimental 
conditions, meaning that techniques developed for polymers are not always directly applicable to liquids. 

This paper gives a description of an experimental method developed for measuring the dynamical expansion 
coefficient of a viscous. As mentioned, the principle is based on the capacitive technique by Bauer et alM^. The 
method is modified in three respects compared to the work of Bauer et ah: 1) The modelling takes into account the 
situation where the capacitor is not full, which is relevant when measuring on a molecular liquid rather than on 
a polymer. 2) The experiment is performed in the time domain using a very fast temperature regulation, which 
gives a dynamical range of more than four decades. 3) The sensitivity is enhanced by using a capacitance bridge 
with a very high resolution. This makes it possible to measure the response following very small temperature steps, 
ensuring that the response is close to perfectly linear. As an application of the technique the paper presents data on 
the glass-forming liquid tetramethyl tetraphenyl trisiloxane (DC704) in the ultraviscous regime. 



In a linear response experiment, the response of a system to an external perturbation is studied. If the perturbation 
is small the output is assumed to be linearly dependent on the input. The formalism to describe this is well known. 
However, different formulations can be used, and the version used in this work when converting the measured time- 
domain response to the frequency-domain response function is maybe not the most common one. The formalism 
used here has the advantage that the time-domain response function and the frequency-domain response function 
have the same dimension and there is no differentiation involved when transforming between the two. This section 
gives a summary of the response function formalism used including a comparison to the standard formalism. 

The fundamental assumption is that the output depends linearly on the input. The most general statement is that 
the change in input dl(t') at time t' leads to a contribution in output dO(t) at time t: 



It is here assumed that the change in output only depends on the time difference (t — t'). Causality implies that 



I. RESPONSE FUNCTIONS WITH CONSISTENT DIMENSIONS 




(1) 



R(t) = Ofor t < 0. 



(2) 



Integrating on both sides of Eq. (T): 




and substituting t" = t — t 



•' and writing l(t) 



dl{t) 



dt 




o 




Changing t" to t': 




(3) 
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If the input is a Heaviside function: 

'<«>-*»<«>-*{! SI^S 



then 

O 



/■OO 

(f) = I / R(t')S(t - t')dt' = I R(t), (4) 
JO 

and it is seen that R(f ) is the output from a Heaviside step input. 

Linear response can also be studied in the frequency domain. In the case of a harmonic oscillating input I(t) = 
j e 1 ( wt +'Pi) r the output O(t) = Ooe l ^ cot+ 't' ^ will be a periodic signal with the same frequency co, but there will be a 
phase shift of the output relative to the input. From Eq. |3j the output is 

Oft) = / R(t , )uol ^ t ^ t ~ t) dt' 
Jo 

= ke'^e't'ito / R{t')e- twt dt' 



= I(t)R(cv), 

where R{oo) is the frequency domain response function, which is given by the Laplace transform of R(t) times ico: 

r°° ■ i 

R(oj) = ico / R(t')e-' wt dt'. 



The linear response relation is often expressed in an alternative formulation where the linearity assumption is ex- 
pressed by 

o(0 = f F (t-t')i(t')dt', 



where pi is sometimes called the memory function, but it is also sometimes called the response function. The use of 
the word response function for ji(t) is somewhat inconvenient because it has a different dimension compared to the 
frequency-domain response function R(to). Substituting again (t" = t — t') and changing t" to t' 



0(t) = / }L {t')l{t-t')dt' 
Jo 



From Eq. @ and |5} we have 



Applying a Heaviside input again 

/• CO 

0{t) = / } i{t')l H{t-t')dt' 
Jo 

= I j*v(t'W. (5) 

R(t)= f F (t')dt', 
Jo 

In the memory function formalism the frequency domain response is again found by inserting a harmonic oscillating 
input. In this case the result becomes 



and therefore 



where the last equality comes from inserting Eq. <(6j> . This expression is formally equivalent to Eq. l|22)l which can be 
shown by integration by parts and by invoking R(t = 0) = 0. However, when converting data in practice Eq. l|22|l 
has the advantages that differentiation of the time domain data is avoided. It is always good to avoid differentiation 
of numerical data because it introduces increased noise. Moreover, if we introduce an "instantaneous" response in 
terms of R(t — > 0) 7^ corresponding to very short times where we can not measures the time dependence of the 
response, then this information would be lost by differentiation. 
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II. PRINCIPLE, DESIGN AND PROCEDURE 

The method requires that there is a simple relation between sample density and dielectric constant. The dielectric 
constant in general has two contributions: atomic polarization and rotational polarization 2 ^. The atomic polarization 
is due to the displacement of the electron cloud upon application of a field. This contribution is governed by the 
microscopic polarizability of the molecule, x (usually called oc, but a is reserved for the expansivity in this paper). 
The atomic polarizability can be assumed to be temperature and density independent in the relevant range. This 
means that the desired simple relation between density and dielectric constant can be obtained when the atomic 
polarization is the only contribution. 

The rotational polarization is due to rotation of the permanent dipoles in the sample. This contribution is relevant 
when the liquid has a permanent dipole moment and mainly at frequencies lower than or comparable to the inverse 
relaxation time of the liquid. The rotational contribution gives the dielectric signal monitored in standard dielec- 
tric spectroscopy. The rotational polarization is temperature-, density- and frequency-dependent, and it is therefore 
non-trivial to relate the density to the dielectric constant when rotational polarization is present. Therefore, in ca- 
pacitative dilatometry it is a contribution one would like to avoid. It is sometimes assumed that the high frequency 
plateau value of the dielectric constant measured in dielectric spectroscopy contains only atomic polarization and 
that it corresponds to the square of refraction index n 2 . However, there is also fast ("glass-like") contribution to the 
rotational part of the polarization. The fast rotational contribution will dominate over the geometric effects even at 
high frequencies if the sample has a high dipole moment. This was demonstrated in Ref. 21. To minimize the rota- 
tional contribution two things are done: 1) Only liquids with very small dipole moment are studied - i.e. liquids in 
which the atomic polarization is dominant at all frequencies and temperatures. 2) These liquids are only studied at 
frequencies much higher than the inverse relaxation time. In the data reported in this paper the measuring frequency 
is 10 kHz and the relaxation time is 100 seconds or more. 

The cell is a capacitor made of circular copper plates of 1 cm diameter and 1 mm thickness, with a 50 f(m spacing. 
The separation is kept by four 0.5 mm x 0.5 mm and 50 /;m thick Kapton spacers. The spacing between the capacitor 
plates is filled with the sample liquid. The thin spacing results in a reasonably large dielectric signal (empty capac- 
itance is 14 pF) despite the small size. The thin spacing moreover makes it possible to heat or cool the sample fast, 
even though the heat diffusion in the sample liquid is slow compared to the heat diffusion in the copper plates. 

The cell is integrated with a microregulator, which is a tiny temperature regulator based on an NTC-thermistor 
(placed in the lower cupper plate of the capacitor-cell), a Peltier-element acting as a local source of heating and 
cooling, and an analog PID-control. The integrated cell and microregulator are placed in our main kryostat. With this 
setup the temperature of the sample can be changed by steps of up to 2 K within less than 10 s and the temperature 
can be kept stable is within a few micro Kelvin over days and weeks. The cell is shown in Fig. 10 and the whole 
system of the main kryostat and the microregulator is described in detail in Ref. [22l 

The principle of the experiment is to make an "instantaneous" step in temperature and subsequently measure 
the capacitance at a fixed frequency as a function of time. From the capacitance we calculate the time-dependent 
expansion coefficient. In order for the temperature step to be "instantaneous" compared to the time scale of the 
relaxation we need the relaxation time to be 100 s or longer. This means that the measurements are performed at or 
below the conventional glass-transition temperature. Nevertheless, it is important to emphasize that, the liquid is 
in equilibrium when the experiment is performed because we wait at least five relaxation times whenever stepping 
to a new temperature before making a measurement. The measurements themselves also must be carried out over 
five relaxation times in order to obtain the relaxation curve all the way to equilibrium. All together, it takes days 
and sometimes even weeks to take a spectrum at a given temperature. This means that the experiment would be 
impossible without the stable temperature control ensured by the microregulator. 

The relaxation time of viscous liquids close to the glass transition is extremely temperature dependent. We there- 
fore need to make small temperature steps in order for the measured response to be linear. This means that the 
change in volume and thereby the measured capacitance is very small, the relative changes in capacitance dC / C are 
of order 10~ 4 . We use an AH2700A Andeen Hagerling ultra-precision capacitance bridge, which measures capac- 
itance with an accuracy of 5 ppm and true resolution of 0.5 attoFarad in the frequency range 50 Hz-20 kHz. The 
capacitance is measured every second at 10 kHz. 

The sample used is liquid at room temperature and the capacitor is filled by letting the liquid suck in using the 
capillary effect. Complete filling is checked by measuring the capacitance before and after filling, comparing to the 
measured dielectric constant measured at the same temperature with a larger capacitor (which is easy to fill). 
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FIG. 2: Schematic drawing of the dielectric measuring cell with the microregulator. The liquid is deposited in the 50 }im gap 
between the discs of the dielectric cell. The Peltier element heats or cools the dielectric cell, depending on the direction of 
the electrical current powering the element. The current is controlled by an analog temperature-control system that receives 
temperature feedback information from an NTC thermistor embedded in one disc of the dielectric cell (reproduced from Ref. |23|) . 



III. GEOMETRY AND BOUNDARY CONDITIONS 



In order to model the relation between the measured change in capacitance and the expansion coefficient some 
assumptions must be made regarding the behavior of the liquid during the experiment. In this section we describe 
these assumption and the arguments on which they are based. 

The capacitor is filled completely at room temperature with a low viscosity molecular liquid. The measuring tem- 
peratures (close to and below the conventional glass-transition temperature) are typically around 100 degrees below 
room temperature for these types of liquids. The cooling makes the liquid contract in the radial direction because the 
distance between the plates is maintained by the spacers (which have a much smaller expansion coefficient). This 
has the consequence that the capacitor is not completely filled at the temperatures where the measurements take 
place. This gives rise to a difference compared to the measurements done on polymers in earlier wor k 13 ' 14 , a differ- 
ence which is taken into account when calculating the relation between the expansion coefficient and the change in 
capacitance in the following section. 

The liquid contracts / expands radially as long as it has low viscosity, but the situation changes when the liquid gets 
ultraviscous. At high viscosities the liquid gets clamped between the plates due to the small distance between them. 
This has the consequence that the liquid can no longer contract/ expand upon cooling /heating by flowing radially, 
but will contract /expand vertically and pull /push the plates changing the distance between them. This effect is 
the basis for the measurement, because the vertical expansion makes the capacitance change, and we calculate the 
expansion from the change in capacitance. 

The distance between the plates is kept by the Kapton spacers at high temperatures (and long times) when the 
sample liquid flows. However, at times where the sample cannot flow, it is the sample, not the Kapton spacers 
which determines the distance. This is true because Kapton has a stiffness 2 ^ of the same order of magnitude as the 
sample (in the GPa-range), but only takes up approximately 1% of the area between the plates. 

The temperature change gives rise to an internal pressure, which is released by pressure diffusion via viscous 
flow. The characteristic time Tj\ QW of the radial flow between two plates of fixed distance I can be estimated by the 
following argument. A temperature step of AT initiates an internal pressure Ap = Kja.„AT in the liquid. This creates 
a radial flow that eventually discharges the surplus volume AV = AToc v nR 1 l. Although the volume flows in the 
radial direction we may as a crude estimation take the volume velocity, V, of planar Pouiseuille flow 2 ^ V = ijp; WZ 3 , 
where L (the dimension in the direction of the flow) can be taken as R, and W (the dimension perpendicular to the 
flow) can be taken as 2nR. The characteristic discharge flowtime then becomes Tjlow = = 6^( j) 2 . The bulk 
modulus and the shear modulus are of the same order of magnitude. It follows that the Maxwell relaxation time is 
roughly given by tm ~ rj/Kj and that Tfi ow tx (R/1) 2 tm- In the experiment we have /=50 }im and R = 5 mm from 
which it follows that the radial flow time is ten thousand times longer than the Maxwell time. The alpha relaxation 
time is roughly given by the Maxwell time, the flow time will be more than ten days when the alpha relaxation time 
is one hundred seconds. This means that the liquid can be considered as radially clamped in the region we study 
(where all relaxation times are longer than 100 seconds). The transition between the radial flow and the clamped 
situation can be seen in dielectric constant when it is measured as a function of temperature, and the observed 
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behavior is consistent with the above estimate. 

The expansion coefficient we study with the boundary conditions described above is not the conventional isobaric 

expansion coefficient, OCp = y IS / because the liquid is clamped in two directions and only free to move in one 

P 

direction. We call this expansion coefficient the longitudinal expansion coefficient, in analogy to the longitudinal 

modulus, (another name for it could be the iso-area expansion coefficient). It is expressed by oc\ = y |K = } Jt / 

where A is the constant area and / is the dimension which is free to respond to the temperature change. The longitu- 
dinal expansion coefficient is related to the isobaric expansion coefficient a„ via the following relation 

1 + 3K T (cv) 

Where G is the shear modulus and Kj is the isothermal bulk modulus, which are both dynamic i.e. frequency or 
time dependent as are the thermal expansion coefficients. 

From this expression we see that oc\ is smaller than Ci„, except at low frequencies (long times, or high temperatures) 
where G = which implies a; = a«. This expression for the longitudinal expansion coefficient is given (but not 
derived) in another equivalent form in terms of the Possoin's Ratio in Refs. 13 17 26 and can be derived from row 3 
of Eq. (53) in Ref . |2/J Also note that there is a total lack of standard notation. Bauer et at use OCp to note the linear 
expansion coefficient, which is the quantity often used to express volume expansion of solids. That is their OCp is 1/3 
of our OCp. The linear expansion coefficient is called cci by Wallace et al2&, while Fukao et alJZ call it a.oo- The quantity 
we call the longitudinal expansion coefficient 0C\ is denoted ocqa (CA for clamped area) by Bauer, «n by Wallace and 
oc n by Fukao (n for normal). 



IV. RELATING THE MEASURED CHANGE IN CAPACITANCE TO a, 
A. Deriving the relation 

In the measurement we perform a small temperature step ST and subsequently measure the capacitance C m as a 
function of time. From the measurements we find the time dependent quantity ^-^p-(t). In the following section 
we show that this quantity is proportional to the expansion coefficient, 0Ci(t), with a proportionality constant that 
depends on €oo and the degree of filling of the capacitor, /, but not on the geometrical capacitance or the distance 
between the plates. 

The starting point is that the only contribution to the high-frequency dielectric constant, £oo, is the atomic polar- 
izability (Sec. [TTJ - We moreover use the Lorentz field™ from which it follows that dielectric constant is given by the 
Clausius-Mossotti relation: 

eco - 1 _ n 
£oo + 2 3eo 

where x is the polarizability of a single molecule, n is the number density of molecules, and eo is the vacuum 
permeability. 

Moreover, we assume that we have a parallel plate capacitor which is partially filled with a dielectric liquid. The 
degree of filling is denoted by / and the measured capacitance is given by 

C m = fe^ + (1 - = [f £oo + (1 - /)] C g (7) 

where Cg = is the geometrical capacitance of the empty capacitor at the given temperature. 
The derivative with respect to temperature is now given by 



_ = + (!-/)]_ +C s f- 



[fe 00 + (l-f)]-J- + Cgf-^. (8) 



Here it is assumed that the liquid does not contract radially at the temperatures (and on the time scale) we consider 

dc i- an H 

~ar anct sr 



(see Sec. [jTTb , thus d//dT = 0. The next step is to calculate and ^f- under the assumption that the area is 
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constant. This was done by Bauer- 3 *^. For completeness we include a detailed derivation as an Appendix. The 
result is 

= -Kie^a, (9) 
where K(eoo) is given by X(eoo) = (eoo - !)(£«, + 2)/3 and 



dT 

Inserting Eq. |(9]l and {TOj) in Eq. ^ yields 

dC m 



dCg 

~ -CgCd . (10) 



dT 



[/e„ + (1 -/)] (-C ?a/ ) - CgfKie^cn 
= -C g [/ eoo + (!-/)+ /K(eoo)] «j • 
Inserting C ? = C m / [/eoo + (1 — /)] and dividing by C m leads to 

1 dC m /eoo + (l-/)+/K(eoo) 



C m dT /eoo + (!-/) 



«2 • (11) 



Isolating finally cl\ gives 



= /goo + (!-/) 1 dc m 

/e co + (l-/)+/K(e o)C m dT 

1 dC„ 



« ; =P(/,eo.)) — -p^, (12) 
Lit! di 



where 



P(/,eoo)- ^ + 



/eoo + (!-/)+ /K(eoo) 



B. The absolute value of a; 

The determination of ol\ and also the uncertainties of the measured value depend on determining correctly the 
proportionality constant P(f, eoo). In order to do so we need to determine the relevant values of / and eoo- To find / 
we use the expansion coefficient and to find the dielectric constant e^ we use the measured empty capacitance along 
with the measured full capacitance. 

The high-temperature expansion coefficient is found 28 to be 0.7*10 -3 K -1 ; at low temperatures we find 29 that it 
is around 0.5*10~ 3 K^ 1 in the long time limit. We use 0.6*10~ 3 K^ 1 as an average value, and find from this that the 
degree of filling is / = 0.95 if the liquid is assumed to contract radially down to 213 K where the relaxation time is 
100 s. The choice of expansion coefficient in the range 0.5-0.7*10 -3 K^ 1 and final temperatures in the range 210-215 K 
makes / change with ±1%. The effect of changing / within this range leads only to ±0.5% changes in P(f, eoo)- 

Isolating the dielectric constant from Eq. 10 gives: 

Cm-CgCl-f) 

= / r K ^ (13) 

From this it is seen that the uncertainty in / also gives an uncertainty in eoo, and this actually has a greater impact on 
the uncertainty of P than the direct effect of the uncertainty on /. Including this effect, the uncertainty in P due to 
uncertain degree on filling is still only ±1%. 

In order to determine eoo from Eq. ((T3l l we need to know the geometric capacitance, Cg. This is found from 
measurements on the empty capacitor at the measuring temperature. We estimate that the uncertainty is ±2% on 
Cg. This estimate is made by comparing measurements made on the capacitor after assembling it anew. The total 
uncertainty on e^ is roughly ±3%, which leads to an uncertainty on P of ±2%. 
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Altogether the uncertainty on P (/,£«>) and therefore on the absolute value of a; is about ±3%. It should be 
emphasized that this uncertainty has no effect on the shape or the time scale of the measured relaxation. This is 
so as long as we stick to linear experiments. For larger temperature steps there will be (at least in principle) some 
second-order effects making P(f, £oo) change during the relaxation because of the change in £co- 

In the modelling of the connection between measured change in capacitance to a; we have not considered the radial 
expansion of the electrode plates. Including this (in the simples possible way) gives rise to an extra additive term 
~C~ T HT m CP " ^ ne s i ze °f this term will be given by the linear expansion coefficient of the electrodes. They are in 
this case made of copper, which at the relevant temperature has a linear expansion of approximately 15 x 10~ 6 K _1 . 
The total measured change in the capacitance is about 50-100 times bigger, thus the effect is small. However, the 
time dependence is different therefore it could be relevant to include this effect in the future. Alternatively we also 
consider shifting to an electrode material with an even smaller expansion coefficient in order to avoid the effect all 
together. 

It should be kept in mind that we have used the Lorentz field. This is an important assumption and the use of an 
other local field, when connecting density with the dielectric constant will change the result. Using the macroscopic 
Maxwell field, will yield the same everywhere, except for X(e TO ) in Eq. I0 which will be given by KMax(e°°) = 
(eoo — 1) instead of the Ki or (eoo) = (eco — l)(£oo + 2)/3. This leads to a 20 % increase in P and the calculated 
numerical value of a/. Again we stress that using another local field will change the absolute values, but will not 
change the time scale or shape of the measured relaxation. 

While none of the above-mentioned things affect the time scale or the spectral shape of the measured relaxation, the 
temperature dependence of £co could in principle affect the temperature dependence of the calculated a;. However, 
this effect is negligible over the 6 degree range studied in the work and P will be considered constant. 

To summarize, the problems discussed in this section can lead to an unknown temperature- and frequency- 
independent scaling of all the measured a;-values. 

C. The shape of the relaxation curve 

In the following we describe the measuring protocol in detail and describe a correction made on the data. We 
moreover use this to give an estimate of the uncertainty on the shape of the relaxation curves reported. 

A main issue is, of course, the first part of the measuring curve where the temperature gets in equilibrium. Fig. Q 
shows details of a single temperature step. It is clearly seen how the target temperature is achieved within less than 
10 s, corresponding to a characteristic time of 2 s. 

Fig. ||4)l a) shows a typical set of temperature steps: a series of up and down jumps are made at the same tempera- 
ture, with variable amplitude. 

Fig. ((Hi b) shows the raw measured capacitance corresponding to the temperature steps in Fig. @ a). Two things 
are worth noticing. First we see the expected rise in capacitance when temperature is decreased. Secondly, we see a 
long time drift of the equilibrium level. At low temperatures where the liquid cannot contract radially it contracts 
vertically. Comparing measurements on the empty capacitor with liquid filled measurements we estimate that the 
expansion coeficient of the liquid is roughly 10 times larger than that of the Kapton spacers. This means that the 
liquid compresses the Kapton. However, on very long times it will be the Kapton which dominates (becuase the 
liquid flows) and the Kapton will therefore slowly relax and press the electrodes apart. We believe that this effect is 
what leads to the long time drift seen in Fig. 1(4)1 b). The drift is subtracted before treating the data, as illustrated in 
Fig. © c) and Fig. ©. 

We make both up jumps and down jumps in temperature and the subtraction of the drift has an opposite effect 
on the two. We can therefore check that the subtraction is made correctly by comparing up jumps and down jumps. 
This is illustrated in Fig. ((6). The superposition of data obtained in up and down jumps also demonstrates that the 
experiment is linear and gives a general estimate of how precise the determination of the curve shape is. 

The comparison of up and down jumps moreover serves to guarantee that the steps are linear. The relaxation time 
is strongly temperature dependent when the liquid is close to the glass transition, and the steps therefore have to be 
very small in order to maintain linear behavior. Smaller steps can be made as well, and the shape of the relaxation 
is maintained, but the curve starts to get noisy because the signal is very small. When we make larger temperature 
steps, we begin to get typical non-linear aging behavior. That is, the relaxation is slower for down jumps than for 
up jumps when the final temperature is the same. The setup is actually well suited for nonlinear experiments also; 
because of the extremely high resolution of the measured quantity we get very well-defined curves and can clearly 
see the nonlinear behavior already at steps of 1 degree. We plan to use the setup for these types of studies, as well, 
but focus in this paper on the linear results. 
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FIG. 3: Zoom on the temperature monitored in the NTC-bead in the lower capacitor plate during the first 40 seconds of at 
temperature step. 




FIG. 4: a) Example of temperature protocol at one temperature. A series of up and down jumps are made at the same reference 
temperature with different amplitude. The temperatures shown are those measured with the NTC-bead in the lower capacitor 
plate. Notice that the smallest jumps are 0.01 K. b) The measured capacitance (blue points). Notice that the relative changes in 
capacitance (dC/ C) for the small jumps is less than 10~ 4 and can still be measured precisely. There is a long time drift in the 
measured capacitance, the dashed line illustrates this background drift and this slope is subtracted of the data before further 
treatment. 

c) The measured capacitance after subtraction of the drift. 



V. DATA 



Fig. 10 shows the expansion coefficient as a function of time at four different temperatures. The data are shown 
for steps made with «0.1 K, except the data at 211 K which are taken with a temperature step of wO.Ol K. This is 
why there is more noise on this dataset. 
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iog 10 (t) [s] 

FIG. 5: The corrected measured capacitance shown on a logarithmic time scale with the temperature change as starting time. The 
temperature steps with the same final temperature (shown in black and blue) all have the same final value of the capacitance. The 
data shown here correspond to the last 8-steps in Fig. {4} 



Fig. |8]l shows all the data from Fig. (0 normalized and superimposed. This illustrates that the measured relax- 
ation obeys time-temperature-superposition (TTS) within the studied (relatively narrow) temperature range. The fit 
shown in figure[7|is a fit to the superimposed curve obtain from the datasets at T=205 K and T=211 K. 

The function used to fit the date is a modified stretched exponential 30 given by: 



oc(t) — Ciao + A# I 1 — exp 



(14) 



In the fit to data we get f> = 0.6 and k = 2.6. The quality of the fit is so good that we have used it as an interpolation 
of the data and used it to calculate the frequency-domain response, which is given by the transformation in Eq. 1(22)1 . 
The transformation is made by making a discrete Fourier transform (using matlabs FFT-procedure) on the fit of the 
normalized curve evaluated in a number of points. The transformed normalized curve is shown in Fig. 10. Here we 
also show an Exponential relaxation which has been transformed using the same algorithm along with the analytical 
Laplace transform. Moreover, the high frequency power law, which corresponds to the exponent of the fit, is also 
shown. 

In Fig. dlfJl l we show the Laplace transformed fit rescaled with amplitudes and time scales in order to show the 
temperature dependence of the frequency dependent thermal expansion coefficient. 



VI. SUMMARY AND OUTLOOK 



We have presented a technique for measuring the dynamical expansion coefficient a. (t) for a glass-forming liquid in 
the ultraviscous range. The experiment is performed on a setup which follows the capacitative principle suggested 
by Bauer et alM: The dynamical range has been extended from 1.5 decade to more than four decades by making 
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FIG. 6: The relative change in C divided by the size of the temperature step. The data shown are the same data as in Fig. ((5) and 
the colors used for each curve is the same. All the curves superpose, which demonstrates that the experiment is linear. Moreover 
it demonstrates that the subtraction of the background drift is successful. 



time-domain experiments, and by making very small and fast temperature steps. The modelling of the experiment 
has moreover been developed. Data is presented on the molecular glass-former tetramethyl tetraphenyl trisiloxane 
(DC704). This data set is to the best of our knowledge the first data on the dynamical expansion coefficient of a 
molecular liquid. 

The technique presented in this paper is based on a principle where the sample is placed in a parallel plate capaci- 
tor such that the sample maintains the spacing between the plates. Thus a change in sample volume in response to a 
temperature change leads to a change of the capacitance. The advantage of this technique is that capacitances can be 
measured with very high precision, and the small density changes associated with linear experiments can therefore 
be determined reliably. One limitation of the technique is that it only works on timescales larger than 10 s, this could 
possibly be overcome by smaller samples and thereby faster temperature control. A more intrinsic limitation is that 
the technique only works for samples with very small dipole moment. For samples with large dipole moment we 
therefore need a complementary technique. 

The measurements of the thermal expansivity is part of a general ambition in the "Glass and Time" group to 
measure different response function of viscous liquids. A unique feature of our techniques is that the all fit into the 
same type of cryostat 22 , ensuring that the absolute temperature of the liquid is the same for all measurements. The 
thermal expansion measurements described in this paper are thus performed in the same cryostat as our shear me- 
chanical spectroscopy 3 -, bulk mechanical spectroscopy 32 , specific heat spectroscopy 3 - 3 - and dielectric spectroscopy 3 -^. 
The properties of liquids close to the glass transition are extremely temperature dependent, and small differences in 
the temperature calibration can lead to rather large differences in the results. Measuring different response functions 
at the exact same conditions therefore makes it possible to analyze new aspects of the viscous slowing down and 
the glass transition. In recent papers we used this to compare time scales of all the different response functions^, to 
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FIG. 8: The data from Fig. l[7} normalized and plotted versus time scaled with the relaxation time (as defined from the fit to the 
modified stretched exponential). The figure demonstrates that the data obey time temperature superposition (TTS). 
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FIG. 9: Illustration of the Laplacetransform. Cyan Diamonds: Laplace Transform of the normalized data (found by brute force 
nummerical integration of the measured points). Blue: Laplacetransform of the normalized fit (see text for details). Red: Expo- 
nential relaxation which has been transformed using the same algorithm as the used for the fit. Black dashed : analytical Laplace 
transform of Exponential relaxation. Blue dashed-dotted line: power law, which corresponds to the exponent of the fit. 

relate linear response to density scaling and to determine the linear Prigogine Defay ratio 3 -^. 
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Appendix I 

The relation between the temperature derivative of the dielectric constant and of the geometrical capacitance with 
the longitudinal expansion coefficient was derived by Baue r 13 ' 14 . For completeness we include a detailed derivation 
in this Appendix. 

The longitudinal expansion coefficient is defined by 



(15) 
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FIG. 10: Laplace transformed fits (see the text for details). The curves are shown a in the dynamical range that roughly corresponds 
to the measurement. 



We start with the temperature derivative of £co, which in this situation is given by: 

dTj {dl) A {dTj A 
so we need and expression for the first term, ^^fp^ • The Clausius-Mossotti relation gives 

xN _ ( eoo - 1) 



(16) 



3e Al (eco + 2) 



(17) 



where N is the total number of molecules, A is the area and / is the thickness such that N/ (Al) is the number density 
of molecules and x the microscopic polarizability of the molecule. 
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We rewrite this to get 

xN 



3e Al 

and take the derivative with respect to I at constant A 



(eco + 2) + 1 



cteoo\ xN ( 8goo \ _, . xN 



which by reinserting Eq. ((T7)l gives 



deoi\ _ (soo — 1) fdeco\ . .l(eoo — 1) 



3/ ; A (eoo + 2) \dl J A y ™ ' ; / ( eco +2)' 
We now isolate ^ ^ in this expression and get 

'9eoo\ 1 (eoo — l)(e«, + 2) 



3/ / 



inserting this in Eq. l fl6l l 

'3eoo\ (e 00 -l)(e 00 +2)l ( dl\ 



dT J T 3 / \dT J A 

which when comparing to the definition of the longitudinal expansion coefficient in Eq. | |T5t can be rewritten as 

'3eoo\ (fioo - l)(eoo + 2) 



dT 



at = -K(eoo)oci 



where the last equality comes from defining X(eco) = ( £ °° !) (£«> +2) 

Now we move on to the temperature derivative of the geometrical capacitance, Cg, which in this situation is given 



The geometrical capacitance itself is given by 



Aeo 
I 



givmg 

'dC g \ Ae 



dl J A ~ P ~ l Cg 

which when inserted in Eq. | |T8t and combined with the definition of the longitudinal expansion coefficient gives 

dC s \ 

3#) =" C ^- (W) 



A 



VIII. APPENDIX II : FD-THEOREM AND THE EXPANSION COEFFICIENT 

This appendix gives the formal definition of the dynamic expansion coefficient, including how it relates to fluc- 
tuations and how the frequency-domain response is related to the measured time-domain response. This is and 
extension of the presentation in Ref. O However, Ref . HH contains a typo as well some definitions which are not 
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precise regarding the absolute levels of the response functions. The precise definitions are important for our use of 
the data in Ref . l36l 

The measured response to an external field, whether in the time domain or in the frequency domain, is directly 
related to the equilibrium thermal fluctuations of the system. This is expressed formally through the fluctuation 
dissipation theorem (FDT), which expressed in the time domain i o 37 < 38 : 

Here R(t) is the response function (see Sec. <JTJ> for a definition) and sharp brackets refers to ensemble averages. A 
is the measured physical quantity (that is the output 0(f) in Sec. (HJl) and B is conjugated to the applied input /field, 
which is called /(f) in Sec. Q. The function (AA(t) AB(O)) is the correlation function, which in the simple case where 
A — B reduces to the auto correlation function. 

Integrating on both sides of Eq. | |20l l and inserting R(t = 0) =0 gives the time-domain response function : 



/„ 



o dt Jo k B Tdt' x v ' K " 



R(t) = ^«AA(0)AB(0))-(AA(f)AB(0))). (21) 

from which it is seen that R(t = 0) = as it should be. The frequency-domain response function, is given by the 
Laplace transform of R(t) times ico: 



R(co) = ico J R(t')e- ia3t dt'. (22) 
Combining this with Eq. (|2T1 l gives the FDT in the frequency domain: 



ico 



R(co) = —p= I (AA(f)AB(O)) — {AA{0)AB(Q))e~ lwt dt 
kbI Jo 

1 — ™ f°° -icot. 



r (AA(0)AB(0))-— J o (AA(t)AB(0))e- Mt dt. (23) 

Considering now a linear experiment where a small temperature step 5T is applied to a system at constant pressure 
at t = 0. Its volume response is subsequently measured as a function of time: 

SV(t) = R(t-t')ST(t'), (24) 

then the response function R(t) is given byR(f) = ^ (see Sector more details on the linear response formalism). 
The time-dependent isobaric expansion coefficient is defined by 

1 SV(t) 

« P (f) " 



V ST 
R(t) 



V 



(25) 



In terms of the FDT (Eq. (HTJ), the relevant fluctuations for ocp(t) are volume and entropy, and the expansion coeffi- 
cient can therefore be expressed in the following way: 

* p (t) = ^((AV(0)AS(0))-(AV(f)AS(0))). (26) 



The frequency-domain response function otp[co) is then (from Eq. | |23l l) 

1 



1 idi r°° 

= ^(AV(O)AS(O)) - yj^f J o (AV(t)AS(0))e-^dt. (27) 
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